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BUCKLING LOADS FOR STIFFENED PANELS SUBJECTED TO COMBINED 
LONGITUDINAL COMPRESSION AND SHEAR LOADINGS} 
RESULTS OBTAINED WITH PASCO, EAL, 

AND STAGS COMPUTER PROGRAMS 


W» Jefferson Stroud, Wlllleni H« Greene, end Melvin Si Anderson 


SUMMARY 


The bucklinq analyses used in PASCO are summarized with emphasis placed on 
the shear buckling analyses, PASCO analyses include the basic VIPASA analysis, 
which is essentially exact for longitudinal and transverse loads, and a smeared 
orthotropic solution that was added in an attempt to alleviate a shortcoming in 
the VIPASA analysis-underestimation of overall shear buckling loads. Buckling 
results are then presented for six stiffened panels loaded by combinations of 
longitudinal compression and shear. The buckling results were obtained with 
the PASCO (VIPASA and smeared orthotropic solutions), EAL, and STAGS computer 
programs. The EAL and STAGS solutions were obtained with a fine finite element 
mesh and, therefore, provide benchmark calculations for the entire range of 
combinations of longitudinal compression and shear loadings considered. 

INTRODUCTION 

Whereas bucklinq analysis procedures that are both fast and accurate have 
been developed for stiffened panels subjected to longitudinal (Nx) and transverse 
(Ny) loadings (for example, VIPASA, refs. 1-3, and BUCLASP, ref. 4), no 
such procedure has been developed for analyzing stiffened panels subjected to 
loadings involving shear (Nxy). VIPASA comes very close to meeting both objec 
tives; however, when the loading Involves shear, VIPASA underestimates the 
bucklinq load for the overall mode-that is. the mode for thich the buckling half- 
wavelength In the direction of the stiffeners is equal to the panel length. 

VIPASA is generally accurate for loadings involving shear when the buckling 


half-wavel«nqth In the direction of the itiffener* Is less than one third the 
panel length. 

Shear buckling analysis procedures In current use Include the following 
modeling approaches; Stiffeners modeled as linked plates with Infinite panel 
length ( VIPASAi ref, 1)5 hinges along plate element connections for local 
buckling and smeared stiffnesses for overall buckling (for example, ref. 5 ); 
approximations using discrete El and GJ stiffeners; and general purpose finite 
difference and finite element approaches (for example, EAL, refs, 6-7, and 
STAGS, refs. 8-9). All of these approaches have shortcomings. The shortcoming 
of the approach used In VIPASA is mentioned In the previous paragraph and Is 
discussed In this report. When stiffnesses are smeared, local deformations that 
contribute to the overall buckling mode are lost. Local deformations are also 
lost when the stiffeners are modeled as El and GJ stiffeners. Finite difference 
and finite element approaches can provide high accuracy by using small meshes; 
however, to obtain accurate results the computation costs are high. 

Because VIPASA Is the buckling analysis In PASCO, an alternate solution 
approach tor predicting overall shear buckling was explored and Incorporated In 
PASCO. That approach, which Is referred to herein as the adjusted analysis. Is 
based on smeared orthotropic stiffnesses. The purpose of this report Is to 
present shear buckling results obtained using PASCO (includes both VIPASA and 
smeared orthotropic solutions), EAL, and STAGS to help evaluate the shear 
buckling analyses In PASCO. These results also provide accurate benchmark 
calculations to evaluate other analysis procedures. 


SYMBOLS 


Values are given In both SI and U.S. Customary Units. The calculations 
were made In U.S. Customary Units. 


E 

El, f.2 

^12 

F 

Fd.O 

•"d.QO 

h,0 

Fs,90 

L 

Nx 


N 

Nw 


xy 


u. V, w 
W 

X, Y. Z 
X, y, z 

X. Y*. h 


X, y^. z^ 


Young's modulus 

Young's modulus of composite material In fiber direction and 
and transverse to fiber direction, respectively 

shear modulus of composite material in coordinate system 
defined by fiber direction 

scalar factor that relates Input loading to eigenvalue (see 
eq. 2) 

value of F for discrete stiffeners with stiffeners In 
x-dl recti on. This Is the standard VIPASA solution 

value of F obtained from equation (3) 

value of F for smeared stiffnesses with stiffeners In 
x-dl recti on 

value of F for smeared stiffnesses with stiffeners In 
y-dl recti on 

panel length 

applied longitudinal compressive loading per unit width of 
panel (see fig. 1) 

applied shear loading per unit width of panel (see fig. 1) 

applied transverse loading per unit width of panel (see 

f1g. 1) 

buckling displacements 

panel width 

overall panel axes 

coordinates In longitudinal, transverse, and lateral 
directions, respectively 

local element axes 

local coordinates for each element making up a panel 


^ strain in x-di recti on 

A buckling half-wavelength 

M Poisson's ratio 

Pl» P2 Poisson's ratios of composite material in coordinate system 

defined by fiber direction, m ■ Mg Ei/Eg 

BUCKLING ANALYSES IN PASCO FOR LOADINGS INVOLVING SHEAR 

Abbreviated descriptions of the buckling analyses used in PASCO are pre- 
sented. Emphasis is placed on the effect of a shear loading. A more complete 
description is presented in references 1 and 11. 

VIPASA Buckling Analysis 

Except for special analysis techniques discussed in a subsequent section 
entitled Adjusted Analysis for Shear, the buckling analysis in PASCO is VIPASA, 
described In references 1-3, 10-12. VIPASA treats an arbitrary assemblage of 
plate elements with each plate element i loaded by N^^ , Ny^ , and N^^ . The 
overall panel coordinate system (X, Y. 2) Is shown in figure 1(a); the local 
plate element coordinate system (X^. Y*. Z*). displacements, loadings, and sign 
convention are shown in figures 1(b) and 1(c). The buckling analysis connects 
the individual plate elements and maintains continuity of the buckle pattern 
across the Intersection of neighboring plate elements. 

The buckling displacement w assumed in VIPASA for each plate element is 

w-fi(y,)cosf 

similar expressions are assumed for the Inplane displacement u and v. The 
functions fj(y^) and allow various boundary conditions to be prescribed 

on the lateral edges of the panel . Boundary conditions cannot be prescribed on 
the ends of the panel . 
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(a) Hat-stiffened panel with applied Nx, Ny, and Nxy loading 
and overall panel coordinate system 


^ X, u 



(h) Local plate element coordinate 
system and displacements. 



(c) Plate element 1 with Inplane 
loading. Directions shown are 
positive for prebuckling loads. 


Figure 1,.- Overall panel coordinate system and local plate element coordinate 
system, displacements, loading, and sign convention. 



For orthotropic plate elementi with no ihear loading, 
solution cos^ provides a series of node lines that are straight, per- 

pendicular to the longitudinal panel axis, and spaced \ apart as shown In 
figure 2. Along each of these node lines, the buckling displacements satisfy 
simple support boundary conditions. For values of \ given by x • L, L/2, 

L/3 L/m, where m Is an Integer, the nodal pattern shown In figure 2 provides 

simple support boundary conditions at the ends of a finite, rectangular panel 
of length L. An example In which X - L/2 Is shown In figure 3. 

Node linns 


"y 

♦ ♦♦♦♦♦♦ 



Figure 2.- Node lines produced by w » fi(yj^) cos^ • 



Figure 3.* Buckling of orthotropic panel under longitudinal loading. 
Mode shown Is x > L/2, 
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For .oUotropic plat. .l«a.nt. anP/or pl.t. .1th a ih.ar loaPtnp. 

f (y ) U not taro. (Bocaui. anltotropy p.n.ra11y hat naplUIhl. for 

lonplav.1.npth hucklinp itod.i anP baeau.o it I. thia. 1onp^.v.l.nc,th laodta 
that ar. troublaaomo, raferanc. to aniaotropy la dropp.d In th. follclng dla- 
cusaion.) Node linot ore skewed and not straight, but the node lines ore 
still spaced X apart as shown In figure 4. Since nod. lines cannot coincide 
with the ends of the rectangular panel, the VIPASA solution for loadings 
involving shear Is accurate only when many buckles form along the panel length, 
in which case boundary conditions at the ends are not Important. An example 
In which X - L/4 Is shown In figure 6. _ 

Ny 

tilt* 



Figure 6v Buckling of panel under shear loading. Mode shown Is X • L/4 


Ai \ approachei L» thi VIPASA buckling analyil* for a panal loadad by 
Njjy can undartatlmati tha buckling load lubitantlally . Ona axplanatlon 1i ai 
follows: As can bt saan In flgura S| tha skawad nodal Hnai g1 van by VIPASA in 
tha casa of shaar do not coincide with tha and edges. Forcing node lines (and, 
therefore, simple support boundary conditions) to coincide with tha anti edges 
produces long-wavelength buckling loads that are, in many cases, appreciably 
higher than those determined by VIPASA. 

In summary, for stiffened panels composed of orthotropic plata elements 
with no shear loading, the VIPASA solution is exact in the sense that it is the 
exact solution of the plate equations satisfying the Kirchoff-Love hypothesis. 
However, for stiffened panels having a shear loading the VIPASA solution can 
be very conservative for the case X ■ L. It Is believed that for this case 
the VIPASA solution provides lower bound. 

Because VIPASA is overly conservative in the case of long-wavelength 
buckling if a shear load is present, an adjusted shear analysis procedure 
based on smeared orthotropic stiffnesses has been incorporated In PASCO. This 
adjusted analysis can be used (at the user's option) for the case x « L. 

Adjusted Analysis for Shear 

In VIPASA, a scalar quantity denoted FACTOR is the unknown In the eigen- 
value analysis. In this report, a quantity Is Introduced that has essentially 
the same meaning as FACTOR. That quantity Is denoted F. The quantities FACTOR 
and r differ only in that whereas FACTOR Is always the solution of an eigenvalue 
analysis in VIPASA and is Identified with the word FACTOR in the VIPASA print- 
out, y may not be the solution of a VIPASA eigenvalue analysis If an adjusted 
shear analysis Is used In PASCO. Otherwise, FACTOR and F are identical. In 


8 


thil rtport, th« ic«1af F 1i d«f1n«d by 


F 



II 


Input 



elqenvilui 


( 2 ) 


Where N^i Ny, and N^y are Inplane loads, A definition of F which Involves 
additional loads Is given In reference 11, 

The rationale presented here for the adjusted analysis approach Is somewhat 
different from that presented In reference 11, The objective of the analysis Is 
to solve the shear buckling problem for the finite panel Illustrated In figure 6, 


For buckling half-wavelength X equal to panel length L, the mathematical model 
solved by VIPASA and the resulting node lines are similar to those Illustrated 
in figure 7. The panel 1n figure 7 Is Infinitely long In the x-directlon. 



Figure 6.- Finite stiffened panel of 
length L and width W, simply 
supported on all four edges, and 
subject to shear load N^y, 


w - - H 



Figure 7.- Node lines given by VIPASA 
for shear buckling with x ■ L, 
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Figure 8.- Node tines for buckling of Infinitely wide stiffened panel. 

It Is assumed that a better approximation to the solution for the finite 
panel would be obtained with the Infinitely wide panel ^ shown In figure 8. 
Unfortunately, the mathematical model Illustrated In figure 8 cannot be analyzed 
with VIPASA because VIPASA requires that the panel be uniform In the direction 
of the Infinite dimension. However, the mathematical model obtained by smearing 
the stiffnesses of the stiffened panel of figure 8 can be analyzed with VIPASA. 
Using the definition of F given In equation (2). let Fg^go be the value of F for 
the smeared orthotropic model rotated 90® as shown In figure 8, 

In an attempt to adjust the smeared orthotropic solution Fg^gg to account 
for the effects of discrete stiffeners, the solution Fg^gg Is multiplied by a 
correction factor F(j^o/*^s,0* Both ^s.O calculatd with VIPASA using 

the Infinitely long model of figure 7. The solution F^^g Is obtained using 
discrete stiffeners - the usual VIPASA solution; Fg,g Is calculated with smeared 

iThat approach was used In reference 13. Stiffeners were treated as El stiffeners. 
10 







Tht corr«ct1on f«ctor It baitd on th« mumptlon thit tht 
rolitlonihip botwnn tho imtirtd solution and tho dlicrtta solution for tht 
Infinittly lonq panti It approxlirtattly tht tamt as that for tht Infinittly widt 
pantl. If that assumption wtrt corrtct, than (Fd,o/Fs,o) Ft, 90 would bt tht 
dlicrtta solution for tht Inflnltnly widt pantl of figurt 8. In rtftrtnct 11, 
that solution is dtnottd Fd,90* 


•'d,90 



%,90 


( 3 ) 


The mathematical models for the four solutions that appear in equation (3) are 
Illustrated In figure 9. 



Figure 9,- Analysis models used to obtain adjusted solution for shear buckling. 


Summary descriptions of the four solution are; 

Fd.O Standard VIPASA solution for A « L. 

Fs.O Infinitely long orthotropic plate solution for the case A ■ L 
obtained with VIPASA by smearing the stiffnesses, 

Fs,90 Infinitely wide orthotropic plate solution obtained with VIPASA by 

smearing the stiffnesses, rotating the direction of the stiffnesses 
by 90®, and interchanging the and Ny loads. The eigenvalue used 
is the lowest of the set for \ - W, W/2, W/3,,,, where W is the 
width of the panel being studied, (fig, 6) 

Fd,90 Obtained from equation (3), Does not Involve an analysis. 
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In PASCO, ir an ad.luittd analysis Is s«lectad for tha A * L modfti the solution 
used for that mode Is the smaller of fd, 9 o and Fg^ciQ. (However, at 

the lower bound for the solution). 

The Input parameter SHEAR is used to Indicate whether the adjusted analysis 
Is to be used for the X > L buckling load. If SHEAR ■ 0, the standard VIPASA 
analysis Is used. If SHEAR # 0, the adjusted analysis Is used and the value of 
the twisting stiffness used In calculating the smeared orthotropic plate buckling 
load 1s the product of SHEAR and the value of the twisting stiffness calculated 
by equations (43) and (44) In reference 11. A value of SHEAR less than 1 Is 
generally appropriate for a panel composed of closed section stiffeners, such 
as a hat-stiffened panel. 

To summarize the various possibilities for F: 

• When the adjusted shear analysis Is used (SHEAR # 0 ar^ X - L), F 1s 
the smaller of F(j^go and F^^go. 

• For all other cases (SHEAR > 0 or X # L), F Is the usual VIPASA 
solution (Fd^o when X ■ L). 

The adjusted analysis Is an engineering approximation, and engineering 
judgment should be used In Us application. For example, the degree to which 
the smeared orthotropic solutions F^^q and F^^go are good approximations to the 
corresponding discrete stiffener solutions depends, In part, on the buckle mode 
shape of the smeared solution. In both cases, the buckle half -wavelength trans- 
verse to the stiffening should be at least 2.S times the stiffener spacing. In 
subsequent sections of this report, calculations for several stiffened panels 
are presented to assess the validity of the adjusted analysis. Based on these 
calculations, recommendations for using the adjusted analysis are presented In 
the section entitled, DISCUSSION OF RESULTS. All of these recommendations 
Involve a reduction In the buckling load predicted by the adjusted analysis. 
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STIFFENED PANEL EXAMPLES 


Six st1ffan«d panals ware analyxed with PASCO and with the general finite 
element structural analysis codes EAL (refSt 6« 7) and STAGS (refs* 8, 9) to 
evaluate the shear buckling analyses (VIPASA and the adjusted analysis) In 
PASCO. Results of these analyses are presented In this section. Four of the 
six panels had blade stiffeners, one panel had hat stiffeners, and one panel 
was a corrugated panel. All panels were 76.2 cm (30 In.) square and had six 
equally-spaced stiffeners. The loadings were combinations of longitudinal 
compression (N^) and shear (N;^)« STAGS results are presented only for the 
pure shear loadings and only for the first four examples— the blade stiffened 
panels. The VIPASA results provide an accurate check of the EAL model for 
pure longitudinal compression, and the STAGS results for the first four examples 
provide an Independent check of the EAL results for pure shear. All standard 
VIPASA solutions (not the smeared orthotropic solutions) Include the effect 
of arilsotropic bending stiffness terms for each plate element making up the 
panel cross section. 

A schematic drawing showing the loading and overall dimensions for the 
six example cases Is shown In figure 10. The manner In which the applied 
loads were distributed over the cross section— the prebuckling stress state— Is 
discussed In reference 11. In particular, the N)( load was distributed assuming 
uniform strain €x of the panel cross section with free transverse expansion 
of each plate element, so that Ny^ was zero. Buckling boundary conditions 
were simple support on all four edges. These boundary conditions are defined 
In figure 10. The panel cross sections were treated as collections of lines 
with no offsets to account for thicknesses. (Offsets are available In PASCO.) 
The first example Is discussed In greater detail than the other examples. 



Bucklinq boiiiuliry conditions «r» simple support on all four etiges 
X - «. I ! u and dw are unrestrained, v - w « ft 
dx 

V ~ ft Wi V and Ow are unrestrained, u - w = 0 

(»y 


Flqure 10.- Loadinq, dimensions, and boundary conditions for stiffened panel 


ExAinpIn 1 - Coiftposite Blad«-St1ffAn«d Panel 


Panel descriptlon t - A repeatinq etement of the compoaite blade-atiffened 
panel Is shown In figure U. Element widths are also shown. The wall con- 
struction for each plate element is given in table I. Only half the laminate 
is defined for each plate element because all laminates are symmetric. Plate 
element numbers ar'e indicated by the circled numbers in figure 11. Fiber 
orientation angles are measured with respect to the X axis, which is parallel 
to the stiffener direction. 

Values of Young's moduli, shear modulus, and Poisson's ratio for the 
graphite-epoxy material used in the calculations for this example are given 
in table II . 



Figure 11.- Repeating element for example 1, composite blade-stiffened panel. 


PASCO input . - Sample PASCO input for this example is shown in figure 12. 
In this input, the loading is “ 1000, ■ 1000 which means that a solution 

is sought for the case * N^y. ICARD input Is included in order to get 
detailed plots of buckling mode shapes. The repeating element shown in figure 
11 was generated with PASCO input. 


IS 


<»**** BXAMPLR 1| COMPOSITE BLADE-8TIPPENED PANBI^ ***** 

$CONDAT 

$ 

$PANBL 

B-2.S, U3S2, 2*5, 

T-.0055, .0055, .0495, .0055, .Oil, 

THET-45, 0, 90, 45, 0, 

KWALUl, 0-1, -1,-1, l.2»3i 
KWALL(l,2)-4,-4,-4,4,5, 

XWAa-1,2,1, 

HCARD-4 ,-4,2, 90,0, 

2,121,4, 

4.5.1.3, -121, 

ICARD-5,1,3,1,-909,0900, 

3. 2. 3.4, 

3 , 3 , 4 , 3 , 

3.4, -909,0900, 

NOBAY-6 , 

ICREP-6 . 

EL- 30, 

MINLAM-30 , 

IBC-1, 

SHEAR-1, 

lP-2, 

NX- 1000., 

NXY-IOOO., 

$ 

SMATER 

E1--19.E6, E2-1.89E6, E12-.93E6, ANUI-.30, RH0-.0571, 

ALFAl ( 1 )— .005E-6 , ALFA2 ( 0-2 1 .8E-6 , 

ALL0W(l,O- 2, .004, -.004, .004, -.004, .01, 

$ 

Figure 12.- Sample PASCO Input for example 1, composite blade-stiffened panel. 


EAL model . - The single finite e'«ment t>pe used In the EAL model for 
this and the other examples Is a four-node, quadrilateral, combined membrane 
and bending element. Both the membrane and bending stiffness matrices for the 
element are based on the assumed stress, hybrid formulation of Plan (refs. 6 and 
14). The buckling or geometric stiffness matrix for the element Is based on a 
conventional displacement formulation that Includes terms allowing Inplane (u 
and V displacements) as well as out-of-plane (w displacements) buckling modes. 
The Plan membrane formulation allows a single element across the depth of a 
blade stiffener to accurately represent Its overall Inplane bending behavior. 
The EAL designation for this element Is E43. 


i 


The finite element grid chosen for the Ul model Is shown In flqure 13. 

Two elements ere used along the depth of the blade, four elements arr used 
between blades, and 36 elements are used along the length, making a total of 
1296 elements and 1369 nodes. A more refined model consisting of three elements 
along the depth of the blade, four elements between blades, and 4B elements 
along the length was used as a check for selected loadings. (In all cases, 
the differences in the results for the two models were negligible). 



(a) Three views of model. 


Figure 13.- EAU finite element model for example 1, composite blade-stiffened 
panel . 
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It 1s recoflnized that the EAL model just described and the STAGS models 
described In the next section are more refined than models used In usual 
enqineerinq calculations. However, since accuracy was an Issue In this study 
and since relatively Inexpensive and accurate procedures do not exist for cal- 
culatlnq bucklinq loads when the loadinq Involves shear. It was decided that 
benchmark calculations that differ from the exact solution by no more than 
approximately one percent were needed. Based on convergence studies and other 
comparisons, It Is believed that the finite element calculations presented In 
this report meet this accuracy requirement. 

STAGS model . - The single finite element type used In the STAGS model of 
this and the other examples Is a six-node, triangular, combined membrane and 
bending element. The element Is based on the Clough-Felippa triangle and has 
a displacement formulation. Midside nodes allow a single element across the 
depth of a blade stiffener to accurately represent Its overall Inplane bending 
behavior. The STAGS designation for this element Is 42?. 

The finite element grid chosen for the STAGS model Is shown In figure 14, 
The total number of triangular elements In the model Is 2000. The total number 
of degrees of freedom In the model Is 21541. A more refined model containing 
four elements between blades was used as a check for selected loadings. 
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Result! . - BuckUnfl rasiiUs obtained with PASCO, EAt, and STAGS fop this 
example are shown In fiflure IS. The curves Indicate PASCO results, and the 
symbols Indicate EAL and STAGS results. PASCO results are discussed first. 

The dotted line at the top of the figure represents the standard VIPASA 
solution for half-wavelength X equal to L/2 and Is discussed subsequently. 

The four complete curves represent the four solution approaches that appear In 
figure 9 and are Identified in the key for figure 15. The solid curve represents 
the standard VIPASA solution for x * L and Is denoted The curve 

below the solid curve and the curve F^^go above the solid curve represent 
orthotropic plate solutions obtained with VIPASA. These orthotropic plate 
solutions are explained In an earlier section entitled Adjusted Analysis for 
Shear. Each point of the curve F^.o represents a solution for x « L where 
L Is the panel length, and each point of the curve Fg^go represents a solution 
for the lowest buckling load of the set X - W, W/2, W/3, ....where W Is the 
panel width. The curve F<j^go represents solutions obtained using equation (3). 

The corners In the Fs^go and F^^go curves that occur at Nx equal to approxi- 
mately 130 kN/m (750 Ib/ln) Indicate a change In mode shape for the Fg^go 
solution. For Nx less than 130 kN/m, the buckling half-wavelength transverse 
to the stiffeners Is equal to 38 cm (15 In.) which Is three times the stiffener 
spacing. For Nx greater than 130 kN/m, the buckling half-wavelength transverse 
to the stiffeners Is equal to 76 cm (30 In) which Is six times the stiffener 
spacing. In both cases, the buckling mode shape meets the requirement that 
for Fs,go to be valid, the buckle half-wavelength transverse to the stiffeners 
must be at least 2,5 times the stiffener spacing. Since the VIPASA solution 
Is exact for loadings Involving only Nx and Ny, both orthotropic solutions Fg^o 
and Fs^go give the same answer for Nxy » 0. For that reason, Fjj^go " Fj^q Tor 
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Figure 15.- Buckling load Interaction obtained with PASCO, EAL, and STAGS for 
example 1, composite blade-stiffened panel. 

As mentioned earlier, the dotted line at the top represents the standard 
VIPASA solution for \ ■ L/2. It is the X • L/2 counterpart of the F(j^o 
solution, which Is for x > L. Buckling modes with shorter wavelengths 
(X B L/3, 1/4,....) have even larger buckling loads. 

For this example, short wavelength modes have higher buckling loads than 
the Fg^go and F^^gg solutions. Therefore, If an adjusted analysis Is pre- 
scribed, the buckling Interaction curve predicted by PASCO Is the Fg^gg 
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curve. (PASCO telectc the lower of F,,,o tod Fo.so for the X - 1. inode, ei 
exp1a1n«d In a previous section.) 

The circular symbols represent EAL results and the triangular symbol at 
Nx - 0 represents a calculation made with STAGS, The STAGS and EAL results 
at Nx ■ 0 agree to within J.S percent. 

For this example, the Infinitely wide orthotropic solution In PASCO 
and the adjusted analysis F^^gg In PASCO give reasonably accurate estimates of 
the solution for all combinations of Nx and N^. For the loading Nx - 0, the 
solution Fs^gg (which, In this case, would be used by PASCO If an adjusted 
analysis were specified) Is about five percent lower than the EAL and STAGS 
solutions. For this same loading, the standard VIPASA solution F^.g Is about 
63 percent lower than the EAL and STAGS solutions. For the loading - 0, 
the standard VIPASA solution and the EAL solution agree to within 0.3 percent. 

Hetalled comparisons and benchmark calculations for six loadings are 
presented In table III. In this table, the quantity denoted FACTOR Is the 
solution In terms of a scale factor for the specified loading. For example, 
for the loading Nx » 350.3 kN/m, Nxy =175.1 kN/m (Nx = 2000 Ib/ln, Nxy - 1000 
Ib/ln) the EAL solution of FACTOR - 0.4764 means that the solution Is 

Nx » 0.4764 X 350.3 » 166.9 kN/m (952.8 Ib/ln), Nxy - 0.4764 x 175.1 - 83.42 kN/m 
(476.4 Ib/ln). 

Finally, the buckling mode shape obtained with EAL for the case Nx ■ 0 Is 
shown In figure 16. The contour plot shown In figure 16(b) shows that the 
buckling half-wavelength transverse to stiffeners Is approximately equal to 
three times the stiffener spacing, which was predicted by the F^.gg solution. 

The buckling mode shapes obtained with PASCO and EAL for the case Nxy - 0 ore 
shown In figure 17. 
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(«) Mode shape obtained with EAL. 


Undefortned shape 

“~— - Buckling mode shape 



(b) Mode shape obtained with PASCO. 


Figure 17,- Buckling mode shape for pure longitudinal compression for example 1 
composite blade-stiffened panel . ^ 


01' ‘■'(v'kf COW't.Pf^f/ 



TABLE I.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 1. COMPOSITE BLADE-STIFFENEO PANEL 


Laytr Numbtri 
Starting With 
Outalda Layer 


Thickneti 


Fiber 

Orientation* 

Deg 


Plate Elements 1 and 3 


Plate Element 2 


.01397 

.01397 

.02794 


TABLE II.- LAMINA PROPERTIES OF GRAPHITE-EPOXY MATERIAL 
USED IN CALCULATIONS 


Symbol 

Value In 
SI Units 

Value In U.S. 
Customary Units 

El 

131. GPa 

19 X 10® psi 

^2 

13.0 GPa 

1.89 X 10® psi 

^12 

6.41 GPa 

.93 X 10® psi 


.36 

.38 





































Exfliiipla 2 - Metal Blade-Stiffened Panel 


Panel description . - A repeating element of the metal blade-stiffened panel 
is shown in figure IH. Element widths and thicknesses are also shown. The 
material properties used in the calculations are E ■ 72.4 RPa (10.5 x 10® psi), 

M .32 . 
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rigure IB,- Repeating element for example 2, metal blade-stiffened panel. 


PASCO input . - Sample PASCO input for this example is shown in figure 19. 

As is the case with all examples in this report, ICARD input is Included to get 
detailed plots of buckling mode shapes. The same numerical results are obtained 
when the modeling Is carried out with HCARD input only. The repeating element 
shown in figure 18 was generated with PASCO input. 

EAL and STAGS models . - The same finite element and grid pattern used in 
the EAL model of the first example are used in this example. The finite element 
and grid pattern used In the STAGS model of the first example are also used in 
this example. 

Results . - The same general approach used for presenting the buckling 
results In the first example is used In this and subsequent examples. Buckling 
results obtained with PASCO. EAL, and STAGS for this example are shown in 
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***** RXAMPtR 2, METAL BLAnE-STlFPRNRD PANEL ***** 

$CQNDAT 

$ 

$PANRL 

R-2.S, 2.'j, 

T-.042, .029, 

THRT- 0, 0, 

KWALL(t,l)*l, 

KWALL(1,2)"2, 

TWALL-I,2,t, 

HCARn-4,-4,2,90,0, 

2,121,4, 

4. 5. 1. 3, -121, 

TCARO"*! , 1 , 3 , 1 ,-909 ,0900 , 

3. 2. 3. 4, 

3, 3, 4, 3, 

3.4, -909,0900, 
lORKP-6 , 

NOHAY-6, 

SHEAR- I , 

KL-30, 

MI NL AM- 30, 

I HO I, 

TP-2, 

NX- toon., 

NLAM-1,2,3, 

S 

$MATKR 

R1-10.SR6, K2-10.3K6, K12-3.9772727E^, ANUI-.32, RHO-.l, 
ALKAl ( I )— .005K-H , ALFA2( I )-2 1 .8K-6 , 

ALI.0W(l,D- 2, .004, -.004, .004, -.004, .01, 

$ 


Mqure 19.- Sample PASCO Input for example ?, metal blade-stiffened panel. 


fitjure <?0. Curves Indicate PASCO results, and symbols Indicate EAL and STAGS 
results. PASCO results are described first. 

A portion of the standard VIPASA solution for \ L/?, shown by the 

dotted line near the top of the flqure. Is below the f^d,90 solutions 

Therefore, If the adjusted analysis approach Is selected, this portion of the 
dotted line Is an upper bound for the PASCO load Interaction curve. 



N,. kNfin 


Figure 20.- Buckling mode Interaction obtained with PASCO, EAL, and STAGS for 

example' 2 . metal blade-stiffened panel. 

The corners In the Fs,go and Fd^go curves that occur at Nx equal to 
approximately 130 kN/m (750 Ib/ln) and 170 kN/m (950 Ib/ln) Indicate a change 
In mode shape for the Fs^go solution. For values of Nx less than about 130 kN/m 
the buckling half-wavelength transverse to the stiffeners Is equal to 25.4 cm 
(10 In.) which Is twice the stiffener spacing. In this case, the buckling mode 
shape does not meet the requirement that for Ps^go valid, the buckling 

ha if -wave length transverse to the stiffeners must be at least 2.5 times the 



stiffener spacing. Fop values of greater than 13Q kN/m but less than 170 
kN/m, the buckling half-wavelength transverse to the stiffeners Is 38 cm (15 
In.) which Is three times the stiffener spacing. For values cf greater than 
170 kN/m, the buckling half-wavelength transverse to the stiffeners Is 76 cm 
(30 In.) which Is six times the stiffener spacing. 

The EAL and STAGS results fall below both the > L/2 line and the 
Fs, 9 o curve, which Indicates that for this case the adjusted analysis approach 
In PASCO Is unconservative. For the « 0 case, an examination of the EAL 
buckling mode shape (fig. 21) shows that the buckling mode Is an overall mode 
(A = L) rather than a A L/2 mode which might have been assumed because 
the A * L/2 solution Is near the EAL and STAGS solutions. One possible 
factor contributing to the error In the smeared orthotropic solution Fg gQ 
near Nx = 0 is. as pointed out In the previous paragraph, that the buckle 



Figure 21.- Shear buckling mode shape obtained with EAL for example 2. metal 
blade-stiffened panel. 
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hAlf-wAvelenqth transverse to the stiffeners Is too short to be a valid solution* 
These results and results for subsequent examples show the danger of depending 
upon a smeared orthotropic solution— oven If the conservative assumption Is 
made that the panel Is Infinitely wide rather than finite. 

Detailed comparisons of solutions from PASCO* EAL* and STAGS for six 
loadings are presented In table IV. The buckle mode shape obtained with PASCO 
for the case » 0 Is shown In figure 22. 


— U ntfef armed shepe 

Bucklln9 mode shape 



Figure 22,- Buckling mode shape for pure longitudinal compression, obtained 
with PASCO for example 2, metal blade-stiffened panel. 


TABLE IV.- BUCKLING LOADS FOR EXAMPLE 2, 
METAL BLADE -STIFFENED PANEL 


Loading 

FACTOR 



iii 

Ibf 

fn 

S* 

Ibf 

In 

^d,0 

%.o 

%.90 

« 

^d,90 

X - L/2 

EAL 

STAGS 

0 

0 

175.1 

1000 

0.3118 

0.3073 

0.9823 

0.9967 

0.8450 

0.8138 

.8179 

70.0 

400 

175.1 

1000 

.2877 

.2838 

.8423 

.8539 

.7742 

.7195 


175.1 

1000 

175.1 

1000 

.2568 

.2536 

.6879 

.6964 

.6849 

.6061 


350.3 

2000 

175.1 

1000 

.2159 

.2137 

.4637 

.4683 


.4444 


875.6 

5000 

175.1 

1000 

.1413 

.1409 

.1975 

.1980 


.1929 


175.1 

1000 

0 

0 

.97103 

.99695 

.97695 

.97103 


.9759 
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EKample 3 - Heavily-Loaded Composite Blade-Stiffened Panel 


Panel description . - A repeating element of the heavily-loaded composite 
blade-stiffened panel is Shown in figure 23, The wall construction for each 
plate element is given in table V. Only half the laminate is defined for each 
plate element because all laminates are symmetric. Plate element numbers are 
Indicated by the circled numbers in figure 23. Values for material properties 
are the same as those used in example 1 and are given in table II. 



Figure 23.- Repeating element for example 3, heavily-loaded composite blade' 
stiffened panel. 


PA.CQ input and EAL and STAGS models . - Sample PASCO input for this 
example is shown in figure 24. The same finite elements and grid patterns 
used in the EAL and STAGS models of the first example are used in this example. 
Results « • Buckling results obtained with PASC0» EAL, and STAGS are shown 
in figure 25. The four complete curves represent the four PASCO solutions for 
overall buckling. The dotted lines indicate standard VIPASA solutions for 
X = L/2, L/4, and L/6. If the adjusted analysis approach Is selected, these 
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***** EXAMPUE 3, COMPOSITE BUDE-STIFPENED PANEL, HEAVILY LOADED ***** 
$CONDAT 
$ 

$PANEL 
CHANGE- 10, 

6-2.5, 1.97, 2.5, 

T-. 00637, .0249, .0416, .00823, .0673, 

THET-45, 0, 90, 45, 0, 

KWALL(l,l)-l,-l,>l,l,2,3, 

KWALL(l,2)-4,-4,-4,4,5, 
lWALL-1,2,1, 

HCARD-4,-4,2,90,0, 

2,121,4, 

4.5.1.3, -121, 

j ICARD-5,1,3,1,-909,0900, 

3 . 2 . 3.4 , 

3, 3.4, 3, 

3.4, -909,0900, 

NOBAY-6 , 

ICREP-6, 

EL- 30, 

MlNLAM-30, 

IBC-1, 

NLAM-1,2,3, 

IP-2, 

NX- 1000., 

$ 

$ MATER 

El-19. E6, E2-1.89E6, E12-.93E6, ANU1-.38, RHO-,0571, 

ALFAl ( 1 )— .005E-6 , ALPA2 ( I )-2 1 .8E-6 , 

ALL0W(1,D- 2, .004, -.004, .004, -.004, ,01, 

$ 

Figure 24.- Sample PASCO Input for example 3, heavily- loaded composite blade- 
stiffened panel . 


dotted lines form an upper bound for the PASCO load Interaction curve. 

For values of less than about 1600 kN/m (9000 Ib/ln) the buckling 
mode shape for the F^^go solution does not meet the requirement that for 
l! Fs, 9 o to be valid, the buckling half-wavelength transverse to the stiffeners 

must be at least 2.5 times the stiffener spacing. For less than 700 kN/m 
(4000 Ib/ln) the buckling half-wavelength transverse to the stiffeners Is 1.5 
times the stiffener spacing. For Nx greater than about 700 kN/m (4000 Ib/ln) 
but less than about 1600 kN/m (9000 lfo/1n) the buckling half-wavelength tranS' 
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Figure 25.- 


Buckllng load Interaction obtained with PASCO* EAL, and STAGS for 
example 3, heavily-loaded composite blade-stiffened panel. 


verse to the stiffeners Is 2.0 times the stiffener spacing. For Nx greater 
than about 1600 kN/m but less than about 1800 kN/m, the buckling half -wavelength 
transverse to the stiffeners Is 3.0 times the stiffener spacing. 

For the loading ■> 0, the smeared orthotropic solutions Fj^o *^*<1 Ps»90 
overestimate the buckling load by about 7.6 percent. For that reason, there 
Is a bulge to the right In the load interaction curve that would be used by 
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PASCO If the edjuited enelyils approach la lel acted, Thli bulge to the right 
repreaenta Incorrect aolutlonas predicted critical valuea of that are 
greater than the value for " 0 ahould be diaregarded. Thia bulging 
phenomenon also exists In other examples In this report. 

The EAL and STAGS results fall below the ^ ■ L/2. L/4, and L/5 curves 
which indicates that, for this case, the adjusted analysis approach In PASCO 
Is unconservative. PASCO errors In the right-hand portion of the Interaction 
curve would be reduced substantially If the bulge discussed In the previous 
paragraph were disregarded. As In example 2, the buckling mode for the EAL 
analysis at Nx ■ 0 Is an overall mode— not a X » L/2 mode. 

Detailed comparisons of solutions from PASCO, EAL, and STAGS for six 
loadings are presented in table VI. The buckle mode shape obtained with EAL 
for the case Nx > 0 Is shown In figure 26. The buckle mode shape obtained 
with PASCO for the case Nxy ■ 0 Is shown In figure 27. 



figure 26.- Contour plot of shear buckling mode obtained with EAL for 
example 3, heavily-loaded composite blade-stiffened panel. 
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Buckling modt ihipi 



Fisure 27.. ’ff! pure longttudinel compreselon, obtained with 

PASCO for example 3, heavily-loaded composite blade-stiffened panel. 


TABLE V.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN EXAMPLE 3 
HEAVILY-LOAOED COMPOSITE BLADE-STIFFENED PANEL * 


Leyer Number. 
Starting With 
Outside Leyer 


Thickness 


cm 


in 


Fiber 

Orientation, 

Deg 


Plate Elements 1 and 3 


0.01618 

0.00637 

45 

.01618 

.00637 

•45 

.01618 

.00637 

-45 

.01618 

.00637 

45 

.06325 

.02490 

0 

.10566 

.04160 

90 


Plate Element 2 


190 

0.00823 

45 

190 

.00823 

-45 

190 

.00823 

-45 

190 

.00823 

45 

45 

.06750 

0 



TABLE VI.- BUCKLING LOADS FOR EXAMPLE 3. HEAVILV- 
LOADED COMPOSITE BLADE-STIFFENED PAnL 


Loading 

FACTOR 


K 


y 

kN 

m 

Ibf 

TfT 

kN 

m 

Ibf 

liT 


m 

^.90 

^d,90 

X • L/2 

EAL 

STAGS 

0 

87.6 

175.1 

350.3 

700.5 

175.1 

0 

500 

1000 

2000 

4000 

1000 

175.1 

175.1 

175.1 

175.1 

175.1 

0 

1000 

1000 

1000 

1000 

1000 

0 

2.9225 

2.6742 

2.4574 

2.0997 

1.5964 

9.9724 

2.8355 

2.6052 

2.4031 

2.0679 

1.5923 

10.7300 

9.2435 
8.0628 
6.7945 
4.8627 
2.6424 
10.7300 
— 

9.5269 

8.2764 

6.9480 

4.9376 

2.6493 

9.9724 

6.6998 

6.0385 

5.4654 

4.5367 

ni 

6.470 
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IxampU 4 « MtUl Bladt-St1ff*n§d Panti With Thin Skin 


Pan>i datcrlotlon . - A rapeatlnq alamant of tho roatal blado-ttlffantd panti 
with thin skin Is shown In figure 28. Element widths and thicknesses are also 
shown, The material properties used In the calculations are E ■ 72.4 8Pa 
(10,6 X 10® ps1)» u ■ .32. Except for the thickness of the skin, this 
panel Is the same as that used in Example 2. 


3.4)4 cm 
<1.332 In.) 


12, 70 cm , 

1 3.0 In.) 








.127 cm 
10.03 In.) 


fi.)3cm 
" (2.5 In.) 






. 0.147 cm 
(0.036 In.) 




Figure 28.- Repeating element for example 4, metal blade-stiffened panel with 
thin skin. 


PASCO input and EAL and STAGS models . - PASCO Input for this example Is 
the same as that given In figure 19 except that T(l) ■ .025. The same finite 
elements and grid patterns used In the EAL and STAGS models of the first example 

are used In this example. 

Results . - Buckling results obtained with PASCO, EAL, and STAGS are shown 
In figure 29, The same general approach used for presenting buckling data In 

previous examples Is used In this example. 

The dotted curve represents the standard VIPASA solutions for X « L/5 and 
L/6. If the adjusted analysis approach Is selected, this dotted curve forms the 
PASCO load interaction curve. If the adjusted analysis approach Is not selected, 
this dotted curve forms the right-hand portion of the PASCO load Interaction 

curve . 
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Flqure 29.* Buckling load Interaction obtained with PASCO, EAL, and STAGS for 
example 4, metal blade-stiffened panel with thin skin* 


For values cf less than about 150 kN/m (850 1b/1n), the buckling mode 
shape for the F^^go solution does not meet the requirement that for F$^go to be 
valid, the buckling half-wavelength transverse to the stiffeners must be at 
least 2.5 times the stiffener spacing. For example, at ■ 0 the buckling 
half -wavelength transverse to the stiffeners Is 1.2 times the stiffener spacing. 
The EAL results are very near the PASCO short-wavelength results Indicated 
by the dotted curve. Except for the case • 0, the EAL results are slightly 
higher than the PASCO short-wavelength results* For the case Nx ■ 0, the lowest 
eigenvalue for both EAL and STAGS appears to be primarily an overall mode 



rathtr thin a local modaf The second eigenvalue from EALt Indicated by the 
square symbol* Is a local fflodei It Is clear from this example that when the 
mode shape for the smeared orthotropic solution Is In gross violation 
of the mode shape requirements (see previous paragraph)* the adjusted analysis 
approach for predicting the overall buckling mode should be used with caution. 

Detailed comparisons of solutions from PASCO* EAL and STAGS for six 
loadings are presented In table VII. 

The buckle mode shape obtained with EAL for the case ■ 0 Is shown In 
figure 30. The buckle mode shapes obtained with EAL and PASCO for the case 
Nxv = 0 are shown 1n figure 31. 



(b) Contour plot. 


FI pure 31 


Figure 30.- Concluded. 


linrtfforniRl ship* 

. Ruckllitq mort» shipf 



(a) Mode shape obtained with PASCO. 


Buckling mode shape for pure longitudinal compression for 
example 4, metal blade-stiffened panel with thin skin. 










Example 5 - Composite Hat-Stiffened Panel 


Panel description . - A repeatinq element of the composite hat-stiffened 
panel Is shown In figure 32, Plate element widths are also shown, The wall 
construction for each plate element Is given In table VMI. Only half the 
laminate 1$ defined for each plate element because all laminates are symmetric. 
Plate element numbers are Indicated by the circled numbers In figure 32, 

Values for material properties are the same as those used In example 1 and 
are g1 ven In table II . 



Figure 32.- Repeating element for example 5, composite hat-stiffened panel. 


PASCO Input and EAL model . - Sample PASCO Input for this example Is shown 
In figure 33. The finite element used In the EAL model Is the same as that used 
In the previous models. The finite element grid chosen for the EAL model Is 
shown In figure 34, There are 36 elements used along the length. One finite 
element Is used for plate elements 1, 3, and 5, and two finite elements are 
used for plate elements 2, 4, and 6. There are 1944 finite elements and 1813 



***** EXAMPLE 5, COMPOSITE HAT-STIPFENED PANEL ***** 
$CONDAT 
$ 

$ PANEL 
GRANGE- 10, 

IBC-1, 

EL-30. 

B-1.3. 2.4, 1.3, 1.3, l.l, 1.3, 

T-. 010315, .009953, .016955, .025383, 

THET-45, 0, 0, 0, 

KWALLd,!)- 1,-1, -1,1, 2, 

KWALL(1,2)- 1,-1, 3, 

KWALLd, 3)- 1,-1, 

KWALLd, 4)- 1,-1, 4, 

IWALL-1,2,1,3,4,3, 

HCARD-4,-16,6,60,0, 

4,-14, 4,-60, 0, 

4.7.14.5.16, 

5.8.1. 2, -7, 3, 

ICARD-5 ,1,2,1 ,-909 , 0900 , 

3. 2. 5. 2, 

3,2.3,14, 

3,3,4, 5, 

3.4.5.16, 

3. 5. 6. 3, 

3,6,-909,0900, 

lCREP-6, NOBAY-6, 

MlNLAM-30 , 

NLAM-1,2, 

SHEAR-. 3, 

IP-2, 

NX- 1000., 

$ 

$MATER 

El-19. E6, E2-1.89E6, E12-.93E6, ANU1-.38, RHO-.0571, 

ALFAl ( 1 )— .005E-6 , ALFA2 (1 )«2 1 .8E-6 , 

ALLOWd,!)- 2, .004, -.004, .004, -.004, .01, 

$ 


Figure 33.- Sample PASCO Input for example 5, composite hat-stiffened panel. 
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Reiulj^. - Buckling results obtained with PASCO and EAL are shown In figure 
35. The same general approach used for presenting buckling data In previous 
examples Is used 1n this example. 

For values of less than about 300 kN/m (1700 lb/1 n), the buckle mode 
shape for the Fg^go solution does not meet the requirement that for F^^go to 
be valid, the buckling half-wavelength transverse to stiffeners must be at 
least 2.5 times the stiffener spacing. In the range cited, the buckling 
ha If -wavelength is twice the stiffener spacing. 
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Figure 35.- Buckling load Interaction obtained with PASCO and EAL for 
example 5, composite hat-stiffened panel. 
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The EAL results (circular symbols) fall below the smeared orthotropic 
solution Fg^ 9 o for the entire ranqe of loadlnps. The EAk result for ■ 0 . 
Is an overall mode. 

Detailed comparisons of solutions from PASCO and EAL for six loadings are 
presented In table IX. 

The buckle mode shape obtained with EAL for the case ■ 0 1s shown In 
figure 36. The buckle mode shapes obtained with EAL and PASCO for the case 
Nyu = 0 are shown In figure 37. 






(a) Oblique view. 


Figure 36.- Shear buckling mode shape obtained with EAL for example 5 
composite hat-stiffened panel* 




(b) Contour plot. 


Figure 37 


Figure 36.- Concluded. 


undeformed shape 

Buckllns mode shape 




(a) Mode shape obtained with PASCO. 


.- Buckling mode shape for pure longitudinal compression for example 5 
composite hat-stiffened panel. 





TABLE VIII.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 5. COMPOSITE HAT-STIFFENED PANEL 


Lf^er Nunteri 
Starting With 
Outside Lever 

Thickness 

Fiber 

Orientation* 

Deg 

cm 

in 

Plate Eletnents 1 and 3 

1 

0.026200 

0.010315 

45 

2 

.026200 

.010315 

-45 

3 

.026200 

.010315 

-45 

4 

.026200 

.010315 

45 

5 

.025281 

.009953 

0 

Plate Element .2 

1 

.026200 

.010315 

45 

2 

.026200 

.010315 

-45 

3 

.043066 

.016955 

0 

Plate Elements 4 and 6 

1 

.026200 

.010315 

45 

2 

.026200 

.010315 

-45 

Plate Element 5 

1 

.026200 

.010315 

45 

2 

.026200 

.010315 

-45 

3 

.064473 

.025383 

0 
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Loading 


0 

300 

600 

1000 

2000 

1000 


xy 


Ibf 

IT 


m 


Ibf 

TfT 


175.1 

175.1 

175.1 

175.1 

175.1 
0 


1000 

1000 

1000 

1000 

1000 

0 


X 


FACTOR 


m 

D 

%.90 

^d,90 

A-L/2 

EAL 

1.2958 

1.1463 

3.5698 

4.0352 

3.6170 

3.192 

1.2085 

1.0796 

3.1933 

3.5744 

3.3815 

2.932 

1.1288 

1.0176 

2.8822 

3.1970 

3.1653 

2.680 

1.0331 

.9419 

2.4229 

2.6576 

2.9046 

2.3268 

.8389 

.7831 

1.4265 

1.5282 

2.3721 

1.4062 

2.9952 

3.0351 

3.0351 

2.9952 


3.0042 























EKtinpIt 6 • CotnpotUa Corrupittd Pinti 


Pw1 datcpiotion i - A ripMtlnq tlpmtnt of tht compoitti coprupatid panal 
It ihown In flgura 30. Elftntnt width! ara alio ihown* The wall conitructlon 
for aach plata alamant 1i pivan In tab! a X. Plata alamant numbari ara Indicatad 
by the circled number# In flgura 38. Only half the laminate 1i defined for aach 
Plata element because all laminates are symmetric. Values for material pro- 
perties are the same as those used In example 1 and are given In table II. 



Figure 38,- Repeating element for example 6, composite corrugated panel. 


PASCO Input and EAL model . - Sample PASCO Input for this example Is shown 
In figure 39. The same finite element used In the EAL model of the first 
example Is used In this example. For the EAL finite element model, shown In 
figure 40. 36 elements are used down the length of the panel, and two elements 
are used across the width of each plate element. 

Results . - Buckling results obtained with PASCO and EAL for this example are 
shown In figure 41. The EAL results fall above the Fj^go curve; therefore, for 
this case, the adjusted analysis approach Is conservative. This conservatism 



***** EXAMPLE 6, COMPOSITE CORRUGATED PANEL 
$CONDAT 

$ 

$PANEL 


IBC*1, 

EL-30, 


B- .73223, 
T-, 005479, 
THET-45, 0, 
KWALL(1,D- 
KWALL(t,2)- 
IWALL- 1,2,1 
HCARD-4,-6,: 


1.46446, 

.016836, 

l.-l,-l,l,2, 

2 , 1 . 

,-450, -1, 


4, “7, 4, 450, "1, 
6,8, 1,6, 3, 7 ,5, 
ICARD-5, 1,2,1,-909,0900, 
3. 2, 3,6, 


3, 3, 4, 3, 


1.46446, 


1.46446, 


3. 4, 5, 7, 

3, 5, 6, 5, 
3,6,-909,0900, 
ICREP-6, NOBAY-6, 
MINLAM-30, 
NLAM-1,2,3,4, 

SHEAR- 1, 

IP-2, 

NX- 1000., 


$ 

$MATER 

El-19. E6. E2-1.89E6, E12-.93E6, ANUl- 

ALPAl ( I )— .005E-6 , ALFA2( 1 )-2 1 .8E-6 
ALL0W(l,D- 2, .004, -.004, .004, -.004, .01, 

9 


>.38, 


***** 


.73223, 


RH0-.0571, 




Flqure 39.- Sample PASCO Input for example 6, composite corrugated panel. 

is true even though the buckle mode shape for the Fj^go solution does not meet 
the mode shape requirements stated earlier. For the case - 0. the buckling 
half-wavelength transverse to the stiffeners is only 1.2 times the stiffener 
spacing. Results from the other examples suggest that a buckle half-wavelength 
of only 1.2 times the stiffener spacing would cause Fj^go to be very unconserva- 
tive. Detailed comparisons of results from PASCO and EAL for six loadings 
are presented in table XI. 
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(b) Oblique view of model. 

Flqure 40,- EAL finite element model for example 6, composite corrugated panel. 


The buckle mode ihepe fop the pure iheer ceie, figure 42, ihowi thit 
buckling Is limited to the area near the lateral edges* It Is possible that 
this mode Is associated with the Inplane boundary conditions (fig, 10) on the 
lateral edges and with the low value of the extenslonal stiffness of the panel 
In the y-dl recti on. The buckle mode shapes obtained with EAL and PASCO for 
the case ■ 0 are shown In figure 43, 


N,. ib/ln. 



N,. hN/m 


Figure 41,- Buckling load interaction obtained with PASCO and EAL for examole 
6, composite corrugated panel. ^ 
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Flqure 42.- Shea'' buckling mode shape obtained with EAL for example 6, composite 
corrugated panel . 


(a) Mode shape obtained with EAL 


UoMarntril ship* 
Biickllhp tunt* ship* 


^ 

^ V \ / \ / V .* V 




(b) Mode shape obtained with PASCO, 


Figure 43*- Buckling mode shape for pure longitudinal compression for example 6 
composite corrugated panel . 


TABLE X.- wall CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 6 I COMPOSITE CORRUGATED PANEL 


Layer Niniber» 
Starting With 
Outside Layer 


Thickness 



cm 

In 

Plate Elements 

li 3i and 5 

0.013917 

0.005479 

.013917 

.005479 

.013917 

.005479 

.013917 

.005479 

.042763 

.016836 



Plate Elements 2 and 4 


.013917 

.005479 

.013917 

.005479 

.013917 

.005479 

.013917 

.005479 




TABLE XI.- BUCKLING LOADS FOR EXAMPLE 6, 
COMPOSITE CORRUGATED PANEL 


0 

87.6 

175.1 

350.3 

875.6 

175.1 


LOADING 







kN Ibf 

m ITT 






.1654 

1.4462 



0011 1.4986 1.248 
8894 1.3077 1.1395 
7780 1.1242 1.0223 
5971 .8352 .7077 

2865 .3703 .2958 

4462 1.4897 1.4918 




























DISCUSSION OF RESULTS 


The beslc conclusion that can be drawn from these calculations Is that a 
buckling solution based on smeared orthotropic stiffnesses should be used only 
with caution. In most of the calculations, the smeared orthotropic solution 
^s,90 overestimated the overall buckling load. Two factors appeared to 
contribute to the error: (1) smeared orthotropic buckling modes with half- 
wavelength transverse to the stiffeners less than 2.5 times the stiffener 
spacing, and {2} shorter wavelength buckling modes, such as X « L/2, which 
fall below or In the vicinity of the Fg^go curve. In all cases, the finite 
element solution for overall buckling falls between the VIPASA solutions for 
X » L and X * L/2. A solution approach for overall shear buckling that 
assumes the buckling mode to be a combination of the first few VIPASA modes 
is being studied. A special procedure is needed to combine these modes in such 
a way that the boundary conditions at the panel ends are satisfied. 

Based on these results, the following rough guidelines for using the 
adjusted analysis are recommended. If the buckling mode shape for the smeared 
orthotropic solution Fj^gg has a half-wavelength transverse to the stiffeners 
that is less than 1.5 times the stiffener spacing, then accept no Increase in 
buckling load above the standard VIPASA solution Fd^g* If Che buckling half- 
wavelength transverse to the stiffeners Is greater than 2.5 times the stiffener 
spacing, then accept no more than 75 percent of the increase above the standard 
VIPASA solution. Make additional adjustments If local modes are present. 

These guidelines are conservative for the six panels considered, but greatly 
underestimate the load carrying ability of the corrugated panel. Finally, for 
structural sizing appllcatons It Is recommended that the adjusted analysis not 
be used. The resulting panels will be light-weight and conservatively designed. 
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CONCLUDING REMARKS 


PASCO' s buckling analysis (VIPASA) 1s rev1ewed« and an Important short- 
coming of that analysis— underestimation of long-wavelength shear buckling 
loads--1s explained. An alternate overall shear buckling analysis based on 
smeared orthotropic stiffnesses Is presented. The alternate analysis Is 
referred to as the adjusted analysis. 

Ruckling calculations were carried out with PASCO. EAL. and STAGS computer 
programs to study the conservatism In the VIPASA solution and to evaluate the 
adjusted analysis. Six stiffened panels were studied: two composite blade- 
stiffened panels, two metal blade-stiff end panels, one composite hat-stiffened 
panel, and one composite corrugated panel. The loadings were combinations of 
longitudinal compression and shear. It Is shown that the adjusted analysis 
may be unconservative. 

The adjusted analysis for overall shear buckling Is an engineering approxi- 
mation and engineering Judgment should be used In Its application. Rough 
guidelines that might be used are as follows. If the buckling mode shape for 
the smeared orthotropic solution Fs,go {part of the adjusted shear analysis) 
has a half-wavelength transverse to the stiffeners that Is less than 1.5 times 
the stiffeners spacing, then accept no increase In buckling load above the 
standard VIPASA solution F(j^ 9 o. If the buckling half-wavelength transverse to 
the stiffeners is greater than 2.5 times the stiffener spacing, then accept no 
more than 75 percent of the Increase above the standard VIPASA solution. It Is 
recommended that the adjusted analysis not be used for structural sizing appli- 
cations. 
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